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Abstract 

We assess Coarse Graining by studying different partitions of the phase 
space of the Baker transformation and the periodic torus automorphisms. It 
turns out that the shape of autocorrelation functions for the Baker transfor- 
mation is more or less reproduced. However, for certain partitions the decay 
rates turn out to be irrelevant, even decay may stop in a finite time. For the 
periodic torus automorphisms, Coarse Graining introduces artificial dumping. 
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1 Introduction 



Coarse graining is a simple way to explain the manifest irreversibility from the 
underlying entropy preserving dynamical laws. The idea which goes back to Ehren- 
fests [0, see also Tolman [Q, is that fine measurements at the fundamental level are 
unattainable, therefore we can and should observe only averages over microstates. 
The averaging is introduced by the observer in addition to the dynamical evolution. 
The resulting loss of information gives rise to entropy increase. In fact all other 
conventional explanations of irreversibility are inventions of extra-dynamical mech- 
anisms to lose dynamical information. Boltzmann's Stossalanzats |3| amounts to 
loss of postcollisional correlations, von Neumann's 0] measurement projection pos- 
tulate amounts to loss of quantum correlations described by the off-diagonal matrix 
elements of the representation of the pure state in terms of a basis of common eigen- 
vectors of a complete system of commuting observables, resulting in the collapse of 
the wave function. Decoherence amounts to loss of information through transfer to 
an unknown and uncontrollable environment ||. 

We shall not discuss here the details of these approaches which involve assump- 
tions additional to the dynamical evolution of a more or less subjectivistic charac- 
ter. For a recent discussion see for example 0. In contradistinction Prigogine and 
coworkers have stressed 0, f| § that irreversibility should be an intrinsic property of 
dynamical systems which admit time asymmetric representations of evolution, with- 
out any loss of information. These time asymmetric representations can be obtained 
by extending the evolution to suitable Rigged Hilbert Space or by intertwining the 
evolution with Markov processes. These intrinsically irreversible systems include 
for example Large Poincare Non-Integrable systems and chaotic systems. We shall 
not discuss further these interesting direction of research. Our objective is to see 
how coarse graining works for simple systems where calculations are controllable 
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and then draw more general conclusions. In Section 2 we define the coarse graining 
projections and the coarse grained evolution. We study coarse graining projection 
of the Baker transformation in Section 3 and of the periodic torus automorphisms 
in Section 4. 

2 Coarse graining projections 

We introduce the basic concepts and notations that will later be used in the de- 
scription of particular systems. Let us consider the configuration space X with the 
measure \x. The partition ( is the finite set {Z k }^ =1 of the cells Z k which satisfy the 
following properties: 

Ufc=i Z k = X, 
fi(Zif]Z k ) = 8ik(JL(Zi), 

where 5i k is the Kroneker symbol. 

Coarse graining is implemented through the averaging projection over the cells 
of the partition £ known also in probability theory | 
tation operator P over (: 

M 

P/(x) = £/ fe lz fc (*) 



I (J as the conditional expec- 

(2) 



Here lz,.( x ) * s ^ ne indicator of the set Z^. 



l, x e z k 
0, x 4l Z k 



and the average value f k of the function f(x) in the cell Z k is 

1 



fk 



dfj,(x)f(x) 



/j,(Z k ) Jz, 

Let us define by Q the orthocomplement of P: 

P + Q = I. 

Q projects onto the fine/detail information eliminated by P. The simplest dynamics 
are cascades of automorphisms S of the phase space X. The observables are usually 



square integrable phase functions in L 2 (X) and they evolve according to the iterated 
Koopman operator 

V : L 2 (X) -> L 2 (X) such that Vf{x) = f(Sx). 

The coarse grained evolution is described by the projected evolution by P @: 

V£ = (PVP) n , n = 0,1,2.... (3) 

In general however 

VP = PVP + QVP. (4) 

Therefore starting with a coarse grained observable Pf, the evolution may regen- 
erate the fine detail QVPf, which is eliminated by the repeated application of the 
projection P. Therefore arbitrary projections P destroy the dynamical evolution. In 
fact as the Koopman operators of chaotic systems are shifts Jl2| , their coarse grained 
projections can be whatever contractive evolution one wishes (Structure theorem of 



bounded operators) |I2], [13| with any desired decay rates. A condition which guar- 
antees the reliability of the coarse grained description is that the evolution V should 
not destroy the cells of the coarse graining partition. In this way minimal dynamical 
information is lost by averaging and the resulting symbolic dynamics is very close 
to the original evolution. Such coarse-graining projections compatible with the dy- 
namics should be distinguished from arbitrary coarse grainings because they are not 
imposed by the external observer but are intrinsic properties of the system. To our 
knowledge only three types of such intrinsic coarse grainings have been proposed. 
Namely, projections onto the i^-partition |7], 14j| , onto the generating partition [|15[ 



of Kolmogorov dynamical systems and onto Markov partitions [16| . 



We shall not discuss further this interesting subject related also to the symbolic 
representation of dynamics [17| but consider arbitrary coarse grainings as so far no 
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a priori reason to decide over the proper natural partition has been proposed to our 
knowledge. 

The influence of the coarse graining to the approach to equilibrium will be studied 
through the decay rates of the autocorrelation functions of the observable phase 
functions /: 

r (n) (f] = Ix(y n f)(x)f(x)dfl(x) 

U) J x P(x)dfi(x) • 1 ) 

The coarse grained autocorrelation function C^p\f) is 

However, in the following we shall see that the autocorrelation functions are not 
sensitive enough to discriminate results of different approaches. In order to analyze 
a system in more detail, we introduce the decay rate at time (stage) n: 

r (n) (f) = - W \U Ul (7) 

This definition is motivated by the following observation. If the autocorrelation 
function decreases (or increases) monotonically at the points n — 1, n, n + 1, it can 
be written as 

C (n) = A(n) + B(n) exp(-/c(n)n) (8) 

at these points. In this case the decay rate is = n(n). As we can usually expect 
that the representation (|8|) is the leading term of the asymptotics for sufficiently large 
n and that A(n), Bin), and n{n) become independent of n for large n, the decay 
rate converges to the decay rate of the system k(oo). Physically speaking, in 
most cases one needs information about the widths and the lifetimes of the system, 
i.e. T^ n \ rather than just the decaying profile C^ n K 



5 



3 Coarse graining the Baker transformation 

We shall apply the coarse graining to the Baker transformation defined [I8| on the 

torus [0, 1] x [0, 1] by the formula 

R , \ / (2x,f), for 0<x< 1/2 

B ^y) = \ (2x-l,y-f), for 1/2<x<1 ■ (9) 

In order to study the applicability of the coarse graining for this transformation, we 

introduce two different partitions ( s and such that fJ>(Z?) = fi(Zj) for all i, j. As 

these partitions have cells of the same measure, the role of geometry is manifested 

in the clearest way. It is worth noticing here that we cannot use for this study an 

one-dimensional map as we cannot choose two different partitions with the same 

measure in one dimension. In Fig. 1, we present both partitions for the number of 

the cells in each direction M = 4. 

The calculation of autocorrelation functions involves integration over each 

cell. The numerical realization of this integration results in loss of accuracy. To 

avoid this problem, we use the fact that the function /w = (P_BP) n / is piece-wise 

constant for n > 1. For such functions, the successive iterations can be written as 

ft +1) ~E fPtiZk fl B{Zi))- (io) 

It is worth noticing that the latter representation is rather effective from the com- 
putational point of view as the sum in Eq. (|T0|) involves very few terms. The only 
remaining integration in Eq. (§) is the calculation of Pf. We use the rectangular 
quadrature formula for this integration. The numerical investigation shows that the 
number of integration points 40 by 40 for each cell is enough to reach the conver- 
gence. 

Now, after the description of the method, we present the results. We start with 
the initial function 

f(x,y) = x + y. 
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With this particular choice, it is possible to calculate the autocorrelation function 
explicitly. A straightforward analytical calculation gives: 

C (n) = ^(l + l n 2- n + ~2-). (11) 

It is interesting to point out that the second term in expression (|TTD corresponds to 
the Jordan block of the second order in the generalized spectral decomposition [19| . 
With expression (pf), the decay rate can be obtained exactly: 

r W=log2-log(l + ^-^). (12) 

The exact expressions ( |TT1 , |T2"|) are used for the comparison with the coarse graining 
results. 

The autocorrelation function and the local decay rates for both partitions with 
M = 200, and analytical results are presented in Fig. 2. The first conclusion we can 
make here is that coarse graining reproduces the correlation function with accuracy 
better than 10 -4 . Hence, if one is interested only in similar integral characteristics, 
coarse graining can be used. 

However, in many cases one needs more detailed characteristics of the evolution, 
for example the decay rates. Here the situation is changed as different partitions 
result in significantly different decay rates. While in this particular case the square 
partition ( s produces results that are surprisingly close to the exact ones, the results 
for the triangle partition differ drastically although they agree up to 11 iterations. 

The decay rates for the partitions with M = 800 are presented in Fig. 3. For these 
partitions the results agree very well in a wider region, up to about 20 iterations. 
However, further iterations again show huge disagreement between the results for 
the triangle partition and the exact results. 

The square partitions may also give irrelevant results. Namely, if we consider 
the rectangle partitions with 2 N subdivisions of the x axis and 2 M subdivisions of 

7 



the y axis, than after N + M — 1 iterations the coarse grained correlation function 
reaches equilibrium exactly. Hence, there is no decay after this number of stages. 
This statement is also true for the square partition when N = M. However, there 
exists rectangle partition with the same cell area that results in the proper decay 
curve. We illustrate this discussion in Fig. 4. Therefore the decay rates are very 
sensitive to the choice of the type of the partition. 

One may expect that despite the big variations of the decay rates, an average of 
them might be stable. To analyze this possibility, we present in Fig. 5 the average 
decay rates 

1 n+10 

An) = J_ (*) 
10 . z l 1 

for n = 10 and n = 20. However also in this case the previous conclusions remain 
valid: while the square partition produces decay rates robust with respect to changes 
of the cell area, the triangle partition gives results which change irregularly with 
respect to changes of the cell area. However, the average decay rate for the 
triangle partition is stabilized when M is rather big. But we cannot restrict ourselves 
to this small number of iterations (n ~ 10) as the decay rates reach their asymptotics 
much later. The results for the ( s show this clearly as the difference of and 
T (20) j g ra th er pronounced. 

We have already mentioned that our results do not depend on the choice of the 
observable function f(x,y). To illustrate this, we present in Fig. 6 the autocorre- 
lation function and the local decay rates for both partitions with M = 800, and 
the initial function f(x,y) = xJx + y. Analytical results are not available for this 
function. One can see that the results have the same qualitative behavior as previ- 
ous ones. Calculations for other initial functions also give the similar results, so our 
conclusions are valid and independent of the initial function. 
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4 Coarse graining the periodic torus automorphisms 

In the previous section we discussed coarse graining applied for the Baker transfor- 
mation. One could see that, despite the problems with the decay rates, the method 
reproduces the autocorrelation functions rather well. Here we shall show that this 
is not always the case, and coarse graining applied to some dynamical systems may 
produce even qualitatively wrong autocorrelation functions. 

Let us consider the periodic automorphisms T of the torus [0, 1] x [0, 1]: 

T ( X \ = ( ( TuX + Tl2 ^ m ° d 1 ^ 
V V J \ ( T 2i% + T 22 y) mod 1 J ' 

The coefficients are integers with the property 

T k = I for some k e N. (13) 

There exist few families of periodic torus automorphisms T with only possible 
periods k — 2, 3, 4, and 6 (Appendix 1). As described in the previous section, we 
construct the autocorrelation functions (|5|f|). Due to the periodicity of the transfor- 
mation T, the exact autocorrelation function is also periodic. However, the coarse 
grained autocorrelation functions may show damping as a result of information loss. 
We present two such examples in Fig. 7 . For the transformation 




we used the partition ( s with M = 20. This type of the partition gives reasonable 
results for the Baker transformation. Here, however, one can see that the auto- 
correlation function approaches equilibrium as if there was damping. Hence, the 
periodicity is not preserved. Moreover, for the simple transformation 
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there exist also partitions giving rise to damping. As an example, we plot in Fig. 7 
the autocorrelation function for the rectangle partition with 5 subdivisions for the 
and 7 subdivisions for the y axis. 



5 Conclusions 

For the Baker transformation, coarse graining reproduces the shape of the auto- 
correlation functions (Figs. 2,6), but decay rates are hardly reproduced (Figs. 2-6). 
Moreover, some partitions like in Fig. 4 give rise to exact approach of equilibrium 
after a finite time while others may give the correct asymptotic behaviour. For the 
periodic torus automorphisms, coarse graining introduces artificial damping of the 
autocorrelation functions. 

The above mentioned conclusions illustrate once more the subjectivistic charac- 
ter of coarse graining. Therefore when employed, coarse graining should be chosen 
with care, for example by selecting partitions intrinsic to the dynamics, as discussed 
at the end of section 2. In fact the answer to the inverse problem of statistical 
physics is that all stationary Markov processes arise as exact projections onto 



the generating partitions of Kolmogorov dynamical systems, in the spirit of Misra- 
Prigigine-Courbage theory of irreversibility 0, [14j| . In fact such intrinsic partitions 
are not defined by any observations but they are objective properties of the dynam- 
ical evolution. 

Concerning the general issue of irreversibility, the possibilities opened by extend- 
ing the dynamical evolution or by intertwining the dynamical evolution with Markov 
processes [0, |8], |9| are a challenging physical and mathematical research direction. 
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Appendix 1. The periodic torus automorphisms. 



The periodic torus automorphisms 



T = 



( 



a b 
c d 



) 



with integer coefficients satisfy the equation 



T k = j 



where k = 2,3,... is the period. As | det T | = 1, we have detT = ±1. Let 
us transform the matrix T to the diagonal or the Jordan form with the unitary 
transformation. As the identity matrix is unchanged, we can conclude that the 
eigenvalues of the matrix T satisfy |Ai )2 | = 1. Let us examine all possible cases 
corresponding to real and complex eigenvalues. 

I. The eigenvalues are real. In this case we have three possibilities: 

1. Ai = A 2 = 1. As the Jordan block cannot be unitarely equivalent to the identity 
matrix, the only possible transformation is the trivial one, T — I. 

2. Ai = A 2 = — 1. As in the first case, the only possible transformation is T = —I 
which has the period k = 2. 

3. Ai = 1, A2 = —1. In this case we have a family of periodic automorphisms. In 
order do describe them, we may use two invariants: the determinant detT = — 1 
and trace trT = 0, and write down the family as 



These automorphisms have the period k = 2. 

II. The eigenvalues are complex. As the matrix has real coefficients, the eigenval- 
ues are conjugated to each other and can be written as Ai = e tlp , A 2 = e~ tip , ip G R. 
Using the invariants of the matrix, detT = 1 and trT = 2 cosy?, we may write the 




with detT = -a 2 - be = -1. 



(14) 
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corresponding family as 

h \ 

with detT = a(2cosy? — a) — be = 1. 



a 6 
c 2 cos if — a 



The cases when <£> = and (p = it give in fact real coefficients and are already 
analyzed. The only remaining cases when the matrix has integer coefficients are: 
1. cosif = 0, ip = 7r/2. Here we have the family 



a b 
c —a 



with det T = -a-bc= 1. (15) 



These automorphisms have the period k — 4. 

2. cos</2 = 1/2, 9? = 7i"/ 3. The corresponding family is 



a b 
c 1 — a 



with det T = a(l — a) — be = 1. (16) 



These automorphisms have the period k = 6. 

3. cos f = —1/2, (p = 27r/3. The corresponding family is 



a b 
c — 1 — a 



with det T = -a(l + a) - 6c = 1. (17) 



These automorphisms have the period k = 3. 

All possible periodic torus automorphisms are given by fll4|) - (]l 
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Figure captions 



Fig. 1. The partitions used in calculations for M — 4. Fig. la shows the square 
partition ( s , and Fig. lb shows the triangle partition £*. 

Fig. 2. The correlation function (the solid line) and decay rates for 
squares (the long-dashed line), triangles (the dot-dashed line), and analytical results 
(the short-dashed line). The number of cells is 200 by 200. 

Fig. 3. Decay rates for squares (the short-dashed line), triangles (the long- 
dashed line), and analytical results (the solid line). The number of cells is 800 by 
800. 

Fig. 4. Decay rates for squares (the dashed line) with number of cells 256 
by 256 (the cell area 1.52588 10~ 5 ), for 254 by 258 rectangles (the dot-dashed line, 
the cell area 1.52597 10 -5 ), and analytical results (the solid line). The second and 
the third lines practically coincide. 

Fig. 5. Average values t$ for n — 10 and n = 20 as the function of number of 
cells M for squares and triangles. The solid line is for for triangles, the long- 
dashed line is for for triangles, the short-dashed line is for r^ ** for squares, and 
the dot-dashed line is for for squares. 

Fig. 6. The correlation function (the solid line) and decay rates for 

squares (the long-dashed line) and triangles (the short-dashed line). The initial 

function is f(x,y) = x^Jx + y. The number of elements is 800 by 800. 

Fig. 7. The correlation function for the periodic torus automorphisms T = 
/ 1 2 \ 

[ ^ with the partition ( s and M = 20 (line 1), and for the T = 

with the partition of 5 intervals for x and 7 intervals for y (line 2). 
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